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Abstract 

Background: In higher plants, inorganic nitrogen is assimilated via the glutamate synthase cycle or GS-GOGAT pathway. 
GOGAT enzyme occurs in two distinct forms that use NADH (NADH-GOGAT) or Fd (Fd-GOGAT) as electron carriers. The goal 
of the present study was to characterize wheat Fc/-GOG/\r genes and to assess the linkage with grain protein content (GPC), 
an important quantitative trait controlled by multiple genes. 

Results: We report the complete genomic sequences of the three homoeologous A, B and D Fd-GOGAT genes from 
hexaploid wheat {Triticum aestivum) and their localization and characterization. The gene is comprised of 33 exons and 32 
introns for all the three homoeologues genes. The three genes show the same exon/intron number and size, with the only 
exception of a series of indels in intronic regions. The partial sequence of the Fd-GOGAT gene located on A genome was 
determined in two durum wheat {Triticum turgidum ssp. durum) cvs Ciccio and Svevo, characterized by different grain 
protein content. Genomic differences allowed the gene mapping in the centromeric region of chromosome 2A. QTL analysis 
was conducted in the Svevo xCiccio RIL mapping population, previously evaluated in 5 different environments. The study 
co-localized the Fd-GOGAT-A gene with the marker GWI\/l-339, identifying a significant major QTL for GPC. 

Conclusions: The wheat Fd-GOGAT genes are highly conserved; both among the three homoeologous hexaploid wheat 
genes and in comparison with other plants. In durum wheat, an association was shown between the Fc/-GOG/\r allele of cv 
Svevo with increasing GPC - potentially useful in breeding programs. 
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Introduction 

Glutamate is a central molecule in amino acid metabolism in 
higher plants. The ot-amino group of glutamate is directly involved 
in both the assimilation and dissimilation of ammonia and is 
transferred to all other amino acids. In addition, both the carbon 
skeleton and ot-amino group form the basis for the synthesis of y- 
aminobutyric acid (GABA), arginine, and proline. Glutamate is 
also the precursor for chlorophyll synthesis in developing leaves 

[!]■ 

As reviewed by Forde and Lea [2], glutamate synthase 
(GOGAT) is the key enzyme involved in the de novo synthesis 
of glutamate. It catalyzes the transfer of the amide group of 
glutamine to 2-oxoglutarate, with the result of two molecules of 
glutamate yielded. 



The history of the discovery of the two enzymes, their structure, 
and gene regulation has been well documented [3,4]. In plants, 
GOGAT enzyme occurs in two forms, depending on the electron 
donor involved in the reaction: it exist as a ferredoxin (Fd) 
dependent (EC 1.4.7.1), and a NADH dependent (EC 1.4.1.14) 
form. 

Both forms are located in plastids, but, while Fd-dependent 
enzyme is usually present in high activities in the chloroplasts of 
photosynthetic tissues, NADH-dependent enzyme is predominant- 
ly located in non-photosynthesizing cells. The role of GOGAT 
enzymes have been well discussed in rice and conifers [5,6]. 
Several studies showed that GOGAT mutations or gene knock- 
outs, with a consequent reduced enzyme activity for both forms, 
seems to be involved in changes in amino acid metabolism 
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[7,8,9,10]. Only a few studies have reported plant gene isolation 
and sequencing, probably due togene lengths and structural 
complexity. For these reasons, the first reported studies on gene 
sequences described the isolation and sequencing of a full-length 
cDNA clone for maize Fd-GOGAT [1 1] and a partial cDNA clone 
for tobacco Fd-GOGAT [12]. The first plant complete genomic 
sequences identified were two genes from Arabidopsis [13]. For the 
Triticeae crops (wheat, barley, rye), no complete sequence has 
been known although partial sequences were reported for barley 
[14] and fragments for wheat [15]. 

Recently, NADH-GOGAT genomic sc-qucnccs has been 
reported for the A and B genomes of tetraploid durum wheat 
{Triticum turgidum ssp. durum) and for the A, B, and D genomes 
of hexaploid wheat (Triticum aeslivum) [16]. Analysis of the gene 
sequences indicates that all wheat NADH-GOGAT genes are 
composed of 22 exons and 2 1 introns. A comparative analysis of 
sequences among di- and mono-cotyledons plants shows both 
regions of high conservation and of divergence. qRT-PCR 
performed with the two durum wheat cvs Svevo and Ciccio 
(characterized by an high and low protein content, respectively) 
indicates different expression levels of the two NADH-G0GAT-3A 
and NADH-G0GAT-3B genes. 

The Fd-GOGAT protein is a monomeric enzyme of 140— 
160 kDa and has been purified from barley leaves as a single 
polypeptide chain containing iron-sulfur and flavin. [17]. Fd- 
GOGAT activity has been mapped to the centromeric region of 
chromosome 2A [18] where we have previously reported a QTL 
for grain protein content (GPC) [19,20]. 

GPC partiaUy determines the nutritional value and the baking 
properties of common wheat (7". aestivum) as well as the pasta- 
making technology characteristics of durum wheat (T. lurgidum 
ssp. durum). GPC is a typical quantitative trait controlled by a 
complex genetic system and influenced by environmental factors 
and management practices (nitrogen and water availability, 
temperature and light intensity). Breeding for high grain protein 
concentration has been one of the main goals for wheat breeders 
for several decades. Simultaneous increases of GPC and grain 
yield has been difficult to achieve because of a negative correlation 
between grain yield and GPC [21,22]. This occurs because grain 
yield is mostly a consequence of starch accumulation, whereas 
protein accounts for less than 10-15% of the grain dr^' weight, and 
any increase in starch accumulation causes a dilution of the 
protein content if it is not accompanied by an equivalent increase 
in N accumulation. So far, any environmental factor that alfects 
grain yield also afiects GPC. 

Recent investigations [23,24,25,26,27,28,29,30,31, ,32,20] indi- 
cated that factors influencing protein concentration in cultivated 
and wild wheats are located on all chromosomes. Heritability 
estimates for GPC ranged from 0.41 [33] to 0.70 [31], depending 
on the genetic material, environment and the computational 
methods. For this reason, the study of this important character 
through traditional methods is difiicult and time consuming. 

In the present work, we d(;termined the DNA sequence of the 
three homoeologous Fd-GOGAT genes in hexaploid wheat, 
analyzed the exon/intron structure, compared the wheat sequenc- 
es to other plants, and studied tetraploid durum grain protein 
content by QTL analysis and identification of candidate genes. In 
particular, we focused our attention on 2A chromosome where the 
Fd-GOGAT gene is located - identifying and characterizing the 
genomic sequence in durum wheat and determining its correlation 
with QTL for grain protein content (GPC). 



Results and Discussion 

Determination of genomic Fd-glutamate synthase 
(GOGAT) gene sequences 

The complete sequences of A, B and D Fd-GOGAT genes of 
hexaploid wheat were obtained by assembling 454 sequences of cv 
Chinese Spring using a partial barley sequence (NCBI accession 
S58774; Gene ID: 548298) as the initial query. The Chinese 
Spring 454 assembly produced one contig comprised of the three 
A, B and D-genome sequences. Three independent contigs were 
then obtained by splitting the three sequences and assigned to A, B 
and D genomes by amplifying genome specific primers [15] on a 
set of nuUi-tetrasomic Unes (NTs) for chromosome group 2 of 
Triticum aestivum cv Chinese Spring [34,35]. Examples of the 
analysis are shown for the A and B genomes in Figure SI in File 
SI. The three homoeologous sequences are given in Text SI in 
File SI. 

Analysis of wheat gene sequences indicates that the wheat Fd- 
GOGAT gene is comprised of 33 exons and 32 introns for all the 
three homoeologous (Figure 1; Figure S2 in File SI). The sequence 
for intron 1 and the part of exon 1 containing the 5-coding portion 
of the exon could cmly be assigned for the D-genome gene 
(Figure 1; Text Si and S2 in File SI). The cause of missing these 
sequences for the A- and B-genome genes is not known, but 
among possible causes are randomness of the shotgun sequencing, 
difliculty in next generation sequencing through the number of 
homopolymers in this region, or that the sequence is so similar in 
all three homoeologous that the software and manual assemblies 
mis-assigned some sequence reads to the D-genome sequence. For 
the A- and B-genome Fd-GOGAT genes, the 5' extension ended 
with the final 1 3 amino acids of the signal peptide (Figure S2 in 
File SI). The 5'-UTR and the remainder of the coding sequence 
for the A-gen()m(- was found in the Transcriptome Shotgun 
Assembly (TRA) data at Genbank (Figure S3 in File SI). 

The exon/intron borders for all three homoeologous wheat Fd- 
GOGAT genes matched the canonical plant borders (GT...AG) 
for all 32 identified introns (Figure 1). Consensus exon/introns 
boundaries were determined using grass EST sequences aligned to 
the genomic sequence and TRA assemblies. The structure of the 
Fd-GOGAT genes is highly conserved among the three homoeo- 
logues. In addition to the same number of exons/introns, no 
differences were observed in exons size and all three encoded 
mature enzymes have identical sequence lengths of 1519 amino 
acids residues (Figure S4 in File SI, Text S2 in File SI). 
Dilferences were observed in the intron sequences - in particular 
differences in length of 1-50 bases in introns 3-5, 8, 14-16, 22-23, 
25-26, 28-31; 50-100 bp in introns 6, 12, 13, 21; and more than 
500 bases in introns 9, 1 1, and 19 (Figure S2 in File SI). The three 
genes vary in length only by intron differences: the Fd-GOGAT-A 
gene is the longest with 1 5,337 bp, the Fd-GOGAT-D is 1 5, 1 76 bp 
long and the gene located on chromosome 2B is the shortest 
atl4,750 bp (all determined from the beginning of the sequence 
encoding the mature protein through the stop codon). The 
s(;quences among the three liomoeologues are highly conserved — 
with only seven amino acid residue differences among the three 
mature polypeptides (Figure S4 in File SI). 

Comparison of Fd-GOGAT genes in otiier species 

Conservation of the Fd-GOGAT sequences was further 
examined by comparing the A-genome wheat Fd-GOGAT 
mature amino acid sequence with a selection of monocot and 
dicot sequences available in public databases (Figure 2). The 
closest match was found with Brachypodium Fd-GOGAT (97%), 
with rice at 94% and maize at 92%. As expected, the dicot 
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Figure 1. Diagrammatic representation of tiie structure of f£/-G06/4 7" genes. The three homoeologous Fd-GOGAT genes of wheat cv 
Chinese Spring are shown. Exons are numbered boxes. Coding sequences are colored, non-coding sequences are uncoloured boxes. Introns are 
intervening lines. Intron 1 is indicated by the dotted line. 
doi:10.1371/journal.pone.0103869.g001 



sequences were more divergent, at 82% for Arahidopsis and 84% 
for poplar. The monocots are all of the same length (1519 amino 
acids), and the dicot examples are of identical length through the 
frrst 1501 residues. The only length different, and the major 
sequence differences between the monocots and dicots is with the 
signal peptide (Figure S3 in File SI) and residues beyond 1501 
(Figure 2). 

As noted above, there are two versions of GOGAT in plants - 
one requiring ferrodoxin and the other NADH. We have 
previously described the wheat NADH-GOGAT genes [16]. 
Figure 3 compares the amino acids sequences of the three wheat 
Fd-GOGAT homoeologous polypeptides with other monocots, 
dicots, and two green algae - the later being the most primitive 
organismal grouping yet found to contain both GOGAT versions. 
As expected, the two higher plant groups form two distinct 
branches for both GOGATs, with the green algae sequences 
distandy related within each of those two branches - consistent 
with the evolutionary distances. Less clarity is found when the 
analysis includes sequences from additional diverse groups; e.g., 
bacteria, cyanobacteria, Archeae, fungi, arthropods (Figure S5 in 
File SI). 

Glutamate synthase gene mapping 

In order to find polymorphisms in the Fd-GOGAT-A gene, the 
"CELl technology" was used in the durum wheat cvs Svevo and 
Ciccio - characterized by high and low grain protein content, 
respectively. 

Allelic variations, SSRs, Indels, and SNPs are the major types of 
molecular markers that can be developed to detect DNA sequence. 
In particular, SNPs are nucleotide variations in the DNA sequence 
of individuals in a population and constitute the most abundant 
molecular markers in the genome. SNPs are widely distributed 
throughout genomes [36], vary in occurrence and distribution 
among species, and are usually more prevalent in the noncoding 
regions of the genome. There are several methods to discover 
SNPs; e.g., Sanger sequencing, DGCE (denaturing gradient 
capillary electrophoresis), denaturating HPLC and enzymatic 
mismatch cleavage. One of the most efficient and reliable 
enzymatic cleavage based method is TILLING (Targeting 
Induced Local Lesion In Genomes) [37,38]. This approach 
requires a treatment of the amplified DNA with CELl endonu- 
clease, or any of a number of single strand endonucleases, after 



heteroduplex formation between the DNA of lines to be investigat- 
ed. CELl is a glycoprotein from celery, and many green plants. It 
cuts a DNA heteroduplex that contains a base-substitution or a 
DNA loop at the 3' most phosphodiester bond of the mismatched 
nucleotides and produces two complementary sized fragments from 
the original amplified product (up to ~ 1,500 bp amplicons). 
Electrophoresis size-separation on polyacrylamide or agarose gels 
is required to visualize any cleaved fragment. 

A set of primer pairs covering the entire Fd-GOGAT-A genomic 
sequence in (Table 1) was analyzed in the two parental lines Svevo 
and Ciccio in order to find SNPs. Out of 19 primer combinations 
amplified and digested with CEL 1 , four of them, (combinations 2, 
4, 5 and 19) showed a digestion pattern in the heteroduplex lane 
but not in the single parental lines digestions, suggesting real 
differences exists between the two cvs (Figure 4). Mutations were 
confirmed by sequencing the fragments (see Materials and 
Method). A partial sequence of Fd-GOGAT-A gene was obtained 
for both cvs Ciccio and Svevo, which allowed finding a total of five 
SNPs and two indels between the two cultivars. Of the five SNPs 
detected, 2 were located in intronic regions (one each for intron 5 
and intron 31) and 3 were located in exonic region (one for each 
exon 6, exon 31, and exon 32). Both indels were located in introns, 
a 2 bp indel in intron 29 and an 8 bp indel in intron 5. 

With the objective to genetically map the wheat Fd-GOGAT-A 
gene, a primer pair (Fd-GF, Fd-GR) was designed in one of the 
genomic region polymorphic between the cvs Ciccio and Svevo 
and analyzes in the RIL population. The linkage map ''Svevo x 
Ciccio" developed by Gadaleta et al. [39] enriched of new DArT 
markers, were used for the gene mapping. The Fd-G marker 
produced, as expected, a single polymorphic fragment of 284 bp 
present in Svevo and lacking in Ciccio. The fragment was also 
physically mapped in the centromeric region of chromosome 2A 
on the bin C-2AS5-0.78 (data not shown). The linkage group 
identifying the 2A chromosome was of 129,9 cM including 31 
markers (12 SSRs, 5 EST-SSRs, and 12DArT) and the GS2-A2 
and Fd-GOGAT-A genes (Figure 4A). 

Relationship between Fd-GOGAT-A and grain protein 
content 

We focused our attention on the Fd-GOGAT gene located on 
chromosome 2A, where several authors found QTLs for GPC not 
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Figure 2. Plant Fd-GOGAT amino acid sequences. Four monocot and two dicot Fd-GOGAT amino acid sequences for the mature protein are 
aligned with Clustal V: Wheat (A genome; present report), Brachypodium (BRADI1G19080), Rice (Os07g46460), Maize (NM_001 1 12223), Arabidopsis 
(CP002688), Soybean. (AK245357). 
doi:1 0.1 371 /journal.pone.01 03869.g002 



linked to pleiotropic efTects of low productivity in difTerent genetic 
materials [18,20,40]. 

The RIL population CiccioxSvevo was evaluated for grain 
protein content (GPC) and grain yield components in five field 
trials in southern Italy. The analysis of variance for GPC revealed 
highly significant differences between the parental lines Ciccio and 
Svevo and among the RIL in each of the field trials, suggesting 
that the RIL population was suitable for studying the putative 
involvement of the GOGA T genes in the control of grain protein 
content. 

QTL analysis reported in Blanco et al. [20] detected ten QTLs 
for grain protein content on chromosome arms IAS, lAL, 2AS 
(two loci), 2BL, 3BS, 4AL, 4BL, 5AL and 6BS. In order to assess 
the putative relation between the Fd-GOGAT-A genes and GPC 
we re-analyzed the RIL data with the Inclusive Composite Interval 
Mapping method [41] in each of the five environments and across 
environments using the "CiccioxSvevo" map data [39] enriched 
with new DArT marker and including the segregation data of the 
new Fd-GOGAT-A marker. 

Among all the putative QTL for GPC in individual environ- 
ments and across environments, only QTLs with LOD&3.0 values 
were considered in the present work. This new QTL analysis 



revealed that the Fd-GOGAT-A gene mapped in the present work 
co-localized with a major QTL for GPC. In particular the Fd- 
GOGAT-A gene co-localised with a GPC-QTL detected on 
chromosome arm 2AS, in the region comprised between the 
markers Xgwm3 72c and the EST-SSR TC8200 1 (including the two 
closer markers Xgwm339 and Xgwm95) significant in two 
environment and across environments. (Table 2, Figures 4 and 5. 
The Svevo allele increasing the trait had a positive additive efiects 
ranging from 0. 13 to 0.27 with a mean value of 0.24. R value, and 
the percentage of phenotypic variance explained by the additive 
effects of the mapped QTL, ranged from 6% to 19.4% between 
environments and the mean was 19.4 across environments. 

Conclusions 

The current report describes the characterization of the bread 
wheat genomic sequence oi Fd-GOGAT genes and their associ- 
ation with grain protein content (GPC). These gene sequences 
were useful to study the grain protein content in durum wheat, a 
quantitative trait controlled by multiple genes and influenced by 
environmental conditions. 
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Figure 3. Phylogenetic tree of GOGAT polypeptides. The mature polypeptides for Fd- and NADH-GOGAT from a selection of plants an green 
algae were aligned with Clustal W and a nearest-neighbour tree generated with MEGA5. All three wheat homoeologues for Fd-GOGAT are in blue and 
red for the NADH-GOGAT version. 
doi:1 0.1 371/journal.pone.01 03869.g003 
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Table 1. Fd-GOGAT-Al specific primer name, sequence and PGR condition used for SNPs detection. 



Primer nsine 


Forward sec|ijence (5'~3') 


Reverse secjuence (5'~3') 


PCR 

conditions 


FdGOGAT/A_l 


CGCCGTCGCTGTTGCCGC 


TGCTGGCCCAATCATTAAACAAGT 


69°C, 


FdG0GAT/A_2 


GCTTGTGGTGTTGGGTTTGTC 


CTCTCTATCAGCTTTCGGCAG 


60°C, 


FdG0GAT/A_3 


CGGTTCCTTTCAATCTATCAGT 


AATCGGATGCTTTAGGGTCAC 


58°C, 


FdG0GAT/A_4_F 


GGAAGCCACAATACAATCTC 


CTAAACAAAAGTAAAGCAGG 


54 C 


DN_FdG0GAT/A_5 


TCTATGAATACTACAAAGGT 


CAACATAAACAAAACCATCT 


54 C, 


DN_FdG0GAT/A_6 


TGACGGAAGGACGGTAGGGG 


GAGGATTGGAAGTTGACAGGCT 


63 C 


DN_FdG0GAT/A_7 


AAACCCTATGGAACTTGGCT 


AAATCGCTGCTTGAAATAAT 


65 °C 


DN_FdG0GAT/A_8 


TCACAAGGGAAGGAGCCAACAT 


TCAGAACGATCGGAAAGCAC 


62°C 


DN_FdG0GAT/A_9 


TACCCTATCAAGTCCTGTCCTG 


CAAACTGATGGGTGCTGAAA 


58°C 


DN_FdGOGAT/A_l 0 


CACGGCCTGCTGTCCCAATAC 


TGGTCACTGTGGGCATCTTG 


68 C 


DN_FdGOGAT/A_l 1 


CATATCTGGCATTGGAAACAT 


AGTTTTCCAGCCTCTTTGCG 


65 °C 


UI\_raLiUL] A 1 / A_ 1 z 


t — rt~f~l — r{"t~ A A Af" hi — TAT/" A 
t_ 1 L]UU 1 LL] AvjAAAL AL 1 A 1 LA 


TTi — rTf~t — r A^ AT^r — n — rhTf~f~ 
1 1 L] 1 1 1 ALjA 1 Liu 1 Li 1 A 1 LiL 


58^C 


DN_FdG0GAT/A_13 


GTCAAAGCTGCTGCACAAAG 


TTAATGGCACTTGTGGCGGT 


60°C 


DN_FdG0GAT/A_14 


CAGATGTTGTTGATGGGTATT 


CCCTTGTGCAATCTTTATCTCA 


58'C 


DN_FdG0GAT/A_15 


TGCATCTGGACGTTTTGGTG 


GCATTTGCCTTAGATACTCCAG 


60°C 


DN_FdG0GAT/A_16 


GGTGTCGGTAAAGCTTGTAG 


TGTGTTTCCGTAAGACCAAG 


58 C 


DN_FdG0GAT/A_17 


CTCAATCAAGCATGCTGGGGG 


CAAGATCAATGTGCTGCGTTT 


62°C 


DN_FdG0GAT/A_18 


TACGAGCCACATTAGCCCAGT 


CACATAATCGTTGGCCTCTC 


60=C 


DN_FdG0GAT/A_19 


GGCAGTCCTTTGGTTGTTTTCT 


TCTTCTTCGCTGGGTGGCA 


62°C 


Each PGR starts with 5 min at 94 C, followed by 35 cycles of 1 min denaturation at 94"C, 1 
elongation at 72"C, and ends with a final elongation of 20 min at 72°C. 


min annealing at the specific annealing 


T ' reported above, and 2 min 



doi:l 0.1 371 /journal.pone.Ol 03869.t001 



The involvement ot Fd-GOGAT-A gene in the control of GPC 
was carried out with three actions: isolation of homoeologous allele 
located on 2A chromosome in two elite durum wheat cultivars 
with different GPC, gene mapping in a segregant population, and 
association between the gene and GPC evaluated in five 
environments. In the present work, we were able to assemble 
the three complete homoeologous genes from the three hexaploid 
wheat genomes using as initial query a partial barley sequence to 
extract and assembly 454 reads from public databases. The three 
homoeologous genes have the same intron/exon structure with 
several differences in both intron and exon. We then used a set of 
aneuploid lines that led us to attribute PCR fragments to the A 
and B genomes. 

In order to quickly screen the gene sequence looking for SNPs 
between our durum cultivars, we followed a "PCR/C£Li 
strategy" similar to a TILLING approach. In this approach, 
CELl cleaves at the site of heteroduplex indicating mismatches in 
the sequences. This allowed us to identify only SNPs between our 
two cultivars avoiding the sequencing of the complete genes. A 
total of five SNPs and two indels were found of which an insertion 
of 8 bp in cv Svevo was used to construct specific markers and 
map the gene in a segregant population (Svevo xCiccio). Mapping 
data of the polymorphic fragment allowed us to identify the locus 
named Fd-GOGAT-A, in the centromeric region of chromosome 
2A. QTL analysis performed with CIM (Composite Interval 
Mapping) confirmed the presence of the marker in a major QTL 
for grain protein content. Several studies localized QTLs for GPC 
on chromosomes of group 2 [18,42,43,44]. The QTL analysis 
carried out in the RIL population CiccioxSvevo, previously 
evaluated for grain protein content in five environments, showed 
that Fd-GOGAT-A co-localized with QTLs for GPC on chromo- 
some arms 2AS - the CIM analysis identified a major QTL with a 



stable effect in two environments and across environments. The 
genomic differences existing between the two cvs might modify the 
final predicted protein functionality or might have a key role in the 
gene regulation and gene expression, determining a different 
GPC. Further investigation are needed to prove this involvement 
through genetic transformation and/or sexual complementation. 

The influence of the group 2 chromosomes on GPC control was 
reported in different genetic material have previously indicated the 
key role of these chromosomes play in the control of GPC. QTLs 
for GPC on 2A and 2B were firstly reported by Joppa and CantreU 
[45] using durum wheat ssp. dicoccoides chromosome substitution 
lines. Prasad et al. [23] reported a protein content QTL on a distal 
region of the chromosome arm 2AS, while Blanco et al. [40] 
identified a GPC-QTL on 2AS near the centromere. More 
recently was found a stable QTL on 2 A {QGPC.USW-A2) 
expressed in three environments and a QTL on 2B [QGpc.usw- 
B3) significant in four of the six environments analyzed [31]. 

The coding sequences of the homoeologous wheat Fd-GOGAT 
genes showed higher conservation among the three than the 
NADH-GOGAT homoeologues previously described [16], with 
only seven amino acid differences among the three Fd-GOGAT 
homoeologues. The conservation extended in comparisons to 
other plant Fd-GOGAT amino acid sequences, with only the C- 
terminal region having major sequence divergence in monocots vs 
dicots. The two forms of GOGAT are clear when the analysis 
compares among plants and green algae (Figure 3). 

The precise role and location of the two forms of GOGAT are 
not completely understood, but the NADH form seems involved in 
development and re-assimilation of ammonia with either the 
cytosol or amyloplasts [46], while the ferredoxin form in a key 
component in photosynthesis and nitrogen fixation within the 
chloroplasts [47]. Both plastid forms are believed to have 
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Figure 4. Genetic map position of Fd-GOGAT. The map positions and QTLs for grain protein content are shown in the side panels. Fd-GOGAT-A 
and GS2-B2 loci, associated with GPC QTLs, respectively on chromosomes 2A and 2B, are highlighted in yellow. Black dots represent centromere. The 
genetic maps of wheat chromosomes 2A and 2B are from the SvevoxCiccio RIL mapping population. 
doi:1 0.1 371/journal.pone.01 03869.g004 



Table 2. Additive gene effects of the detected QTL for grain protein content in the 2AS region (flanked by the markers TC82001 
and Xgwm372c and including FdCOCAT-Al) in the SvevoxCiccio RIL population grown at five environments. 



Environment 


Effect' 


LOD 




Valenzano 2006 


0,13 


1,6 


6,4 


Foggia 2006 


0,27 


4,7 


17,4 


Gaudiano 2006 


0,27 


3,5 


13,3 


Valenzano2007 


0,18 


2,8 


10,7 


Foggia 2007 


0,20 


1,5 


6,0 


Across environments 


0,24 


5,3 


19,4 



^Effect: positive additive effects are associated w/ith an increased effect from Svevo allele. 
R^: Percentage of phenotypic variance explained by the additive effects of the mapped QTL. 
doi:l 0.1 371/journal.pone.Ol 03869.t002 
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Figure 5. Grain protein QTLs. LOD score scan on chromosome 2A for QTLs associated with grain protein content. The significant scan for QTLs for 
each environment: A) mean across environments; B) Foggia 2006; C) Gaudiano 2006. The position and the name of molecular markers are shown on 
the chromosome along the horizontal axis. The LOD score scan was obtained by ICIM with highlights the markers used as cofactors. LOD stands for 
logarithm of the odds (to the base 10). A LOD score of three or more is generally considered significant - a LOD score of three means the odds are a 
thousand to one in favour of genetic linkage. 
doi:1 0.1 371/journal.pone.01 03869.g005 



originated from endosymbiotic cyanobacteria [48]. The phyloge- 
netic tree of Figure S5 in File S 1 initially suggests cyanobacteria 
{Cyanohacterium, Cyanoihece, Leplolynghya) may possess both 
forms, but closer examination of the sequences finds different 
accessions of each of those three species contains only a single gene 
which can cluster with either the Fd or NADH GOGAT forms 
(not shown). When combined with the anomalies in annotations 
and the complex branching in Figure S5 in File SI, a more 
detailed analysis of GOGAT genes in all phyla is needed to better 
understand GOGAT gene evolution. 

The comparison of a set of plant Fd-GOGAT genes (wheat, 
Brachypodium, rice and maize) suggests regions of greater 
sequence and structure conservation likely related to critical 
enzymatic functions and metabolic control. The higher identity 
was ()bs(T\-c(i betwc-en wheat and Brachypodium serjuenccs, as 
expected due to the genetic closest between the two species. 

Although the two forms of GOGAT, Fd and NADH, catalyze 
the same reaction, the gene structures and their roles in plant 
metabolism are not identical. The two forms have detectable 
conservation in amino acid sequence up to the point where the 
NADH-GOGAT form encodes a pyridine nucleotide-disulfide 
oxidoreductase domain at the C-terminus (arrowhead in Figure 2; 
dot blot in Figure S6 in File SI) - indicating additional metabolic 
aspect of the NADH-GOGAT gene. 

In higher plants, ammonium, whether resulting from nitrate 
assimilation or from secondary sources, is first incorporated into 
glutamine in a reaction catalysed by glutamine synthetase, and 
then glutamate synthase (glutamine:oxoglutarate amidotransfer- 
ase; GOGAT) catalyses the combination of glutamine with 
2:oxoglutarate to form two molecules of glutamate, one of which 
serves as substrate for GS, while the other one is available for 
transport, storage or further metabolism. These two reactions form 
a cycle referred to as the GS/GOGAT pathway [49]. Tissues and 
subcellular localization of both GS and GOGA T genes, as well as 
their different expression level during plant growth, resulting in 
different enzymatic activity along all phenological stages, deter- 
mine a specific synergy between the two genes [47] . In particular 
cytosolic GSl is involved in the pathway with NADH-GOGAT, 
while plastidic GS2 works preferentially with Fd-GOGAT. This 
seems to be confirmed in our genetic material, where both genes 
Fd-GOGAT-A and GS2 were associated with QTLs for GPC. In 
fact, a GS2 gene on chromosome 2B (GS2-B2) was mapped and 
found to be involved in GPC control (Figure 4b) [19]. So, from 
these studies, it's suggested that Fd-GOGAT-A gene works in 
synergy with GS2-B2 gene. This hypothesis could be confirmed by 
further analysis such as gene expression analysis. 

The present work determined the DNA sequences of the three 
Fd-GOGAT genes in hexaploid wheat and identified an allele for 
the increment of GPC in durum wheat and open the way to further 
investigation using the forward and reverse genetic approaches that 
have been successfully used to validate the role of Fd-GOGAT-A 
genes for grain production both in rice and maize [43] . 

Materials and Methods 

Plant material 

Two durum wheat cultivars (Ciccio and Svevo) were used to 
investigate the relation of Fd-GOGAT genomic sequences to GPC. 
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The durum cultivars are parents of a mapping population 
represented by a set of 120 recombinant inbred lines (RILs) 
[.39]. The two parents were chosen for differences in important 
qualitative and quantitative traits; i.e., grain yield components, 
grain protein content, yellow pigment, and adaptive traits. 

Genomic DNA was isolated from fresh leaves using the method 
described by Sharp et al. [50] and subsequendy purified by 
phenol-chloroform extraction. DNA amplifications were carried 
out as described by Gadaleta et al. [39]. 

NuUi-tetrasomic lines (NTs) of Triticum aestivum cv Chinese 
Spring [34,35] were used to physically localize Fd-GOGAT 
markers to chromosomes. Chinese Spring di-telosomic lines [51] 
were used for the assignment of markers to each chromosomal 
arm. Physical locaticm on chromosome bins of each PGR fragment 
was obtained using a set of wheat deletions lines dividing the A and 
B genome chromosomes in bins (kindly provided by B. S. GiU, 
USDA-ARS, Kansas State University) [52]. 

Genes in cv CPiinese spring 

To isolate the complete sequences of Fd-GOGAT genes in bread 
wheat, we used the cDNA sequence of a partial barley Fd- 
GOGAT mRNA (Gene S58774) [14]. This sequence was used as 
the initial query probe to extract matching 454 genomic sequences 
of cv Chinese Spring publicly available (http://www.cerealsdb.uk. 
net/ CerealsDB/Documents/DOC_search_reads.php), which 
were then assembled using the SeqMan module of the Lasergene 
software (DNAStar, Inc.). 

Fd-GOGAT sequences in the Italian durum wheat cvs 
Ciccio and Svevo 

By using Ohgo Explorer and PrimerS (http://frodo.wi.mit.edu/ 
primerS/) software, a set of genome specific primer pairs were 
designed for the distinct Fd-GOGAT hexaploid sequences 
previously obtained: Chinese Spring 2A and 2B Fd-GOGAT. A- 
genome specific primer pairs were used to amplify target DNA 
from both Ciccio and Svevo parental lines using PGR condition 
reported in Table 1 . 

Single PCR fragments were directly purified with EuroGold 
Cycle Pure Kit and sequenced. Multiple PCR products were first 
cloned into the pCR4-TOPO vector (Invitrogcn, Cloning Kit) 
following the manufacturer's instructions, and then each fragment 
was sequenced (http://www.bmr-genomics.it/). 

Sequences assembly were carried out using CodonCode Aligner 
and Geneious software. Gaps and uncertain sequence were 
resolved by primer walking. Regions of less coverage or ambiguous 
reads were rechecked with primers designed to cover those 
regions. 

Digestion with CELT and revelation fragments 

In order to discover and map unknown mutations between the 
genomic sequences of the two varieties Ciccio and Svevo, Single 
Nucleotide Polymorphisms (SNPs) were detected using the 
Surveyor nuclease kit (Transgenomic, Inc.), following manufac- 
ture's instruction. This approach requires a treatment of the 
amplified DNA with CFLI endonuclease, or any of a number of 
single strand endonucleases, after heteroduplex formation between 
the lines to be investigated. Surveyor nuclease cleaves with high 
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specificity at tlie 3' side of any mismatch site in both DNA strands, 
including all base substitutions and insertion/deletions up to at 
least 12 nucleotides. 

Heteroduplex formation, CEL I digestion and gel analysis 

An aliquot of each PGR product was used for hybridization to 
form heteroduplexes between the parental lines, following thermal 
cycler program: 95°C for 2 min; loop 1 for 8 cyties (94°C for 20 s, 
73°C for 30 s, reduce temperature 1°C per cycle, ramp to 72°C at 
0.5°C/s, 72°C for 1 min); loop 2 for 45 cycles (94°C for 20 s, 65°C 
for 30 s, ramp to 72°C at 0.5°C/s, 72°C for 1 min); 72°C for 
5 min; 99°C for 10 min; loop 3 for 70 cycles (70°C for 20 s, 
reduce temperature 0.3°C per cycle); hold at 8°C. After annealing, 
DNA has been treated with Surveyor nuclease to cleave 
heteroduplexes by adding 0.2 |J.l of enzyme and 1.3 |a,l of Buffer 
IX. The digestion step was done at 5°C for 90 minutes and 
stopped immediately by adding 5 |j.l of 0.225 EDTA and 2 \xl of 
bromophenol blue loading dye per sample and mixing thoroughly. 

PGR fragments of over 1000 bp were analyzed on 3% 
polyacrylamide gels. The 3% polyacrylamide gel was made with: 
15 ml 5xTBE, 120 1 dH20, 11 ml 40% bis-acrylamide, 110 |il 
TEMED and 1 ml 10% APS. We used 100-lane \'ertical 
electrophoresis syst(-m (CBS Sci(;ntific, D("l Mar, CA, USA). Th(^ 
images were analyzed manually on PowerPoint (Microsoft Corp., 
Seatde, WA, USA). 

In order to confirm the polymorphisms within genome specific 
genes, the heteroduplex hybridization digestion pattern was 
compared to the ones obtained in each parental lines. Moreover, 
the PCR product giving a digestion pattern after CELl treatment 
were re-amplified and sequenced (http://www.bmr-genomics.it/). 

Development of Fd-GOGAT specific markers, mapping 
and correlation with grain protein content 

The Fd-GOGAT serjuences of the two cvs Svevo and Ciccio 
were aligned using ClustalW from EBI website to identify 
polymorphisms. The marker Fd-G (Forward 5'-GGAAAAGAAG- 
CAGGGCAGATA-3', Reverse 5'-TAGGTCCGTTGGGCAA- 
TACAT-3') for Fd-GOGAT-A, was designed with Oligo Explorer 
software in the polymorphic regions. The polymorphic marker was 
mapped in the "Svevo xCiccio" mapping population. The 
observed segregation ratio for the marker was tested by chi- 
square analysis for deviation from the expected 1:1 ratio. The 
Unkage analysis was performed by JoinMap v. 4.0 [53] and the 
Kosambi mapping function was used to calculate map distances 
[54]. 

Grain protein content (GPC) and yield components were 
evaluated in the RIL population "Svevo x Ciccio" in five different 
environments (Valenzano 2006, Gaudiano 2006, Foggia 2006, 
Valenzano 2007, Foggia 2007). QTL analysis was performed 
following the procedure indicated by Blanco et al. [20] . 

DNA sequence analysis 

DNA sequences were analyzed using the Seqman and Megalign 
modules of the Lasergene software (DNAstar, Inc.), and MEGA5 
[55]. 

References 

1. Yaronskaya E, Vershilovskaya I, Poers Y, Alawady AE, .'\\crina N, ct al. (2006) 
Cytokinin effects on tetrapyrrole biosynthesis and photosyndictic activity in 
barley seedlings. Planta 224: 700-709. 

2. Forde BG, Lea PJ (2007) Glutamate in plants: metabolism, regulation, and 
signalling. J Exp Bot ,58: 2.=i39-23,'i8. 

3. Lea PJ, Miflin BJ (2003) Glutamate synthase and the synthesis of glutamate in 
plants. Plant Physiol Bioch 41: 555-564. 



Supporting Information 

File SI Supporting figures and text. Figure SI. Chro- 
mosome mapping of the Fd-GOGAT genes on chromo- 
somes 2A and 2B. Genome specific markers were amplified 
from durum cultivars and hexaploid wheat cv Chinese Spring 
genetic stocks. A) The A-genome specific marker amplified in cvs 
Svevo, Ciccio, and Chinese Spring and sets of Chinese Spring 
nuUi-tetrasomic deletion lines for chromosome group 2. The 
350 bp fragment was absent in the nulli-2A-tetra-2B line, as 
indicated by arrows and confirming the localization on chromo- 
some 2A. B) The B-genome specific marker amplified in cvs, 
Ciccio, Svevo, Chinese Spring and sets of nullitetrasomic deletion 
lines for chromosome group 2. The 450 bp fragment was absent in 
the nulli-2B-tetra-2A Une, as indicated by arrows and confirming 
the localization on chromosome 2B. Figure S2. Alignment of 
wheat FD-GOGAT genes. Alignment of wheat cv Chinese 
Spring A, B, and D genome Fd-GOGAT genes from the beginning 
of exon 2 through the stop codons. Exons are indicated by red 
brackets and exon number above the aligrmients. The blue line 
indicates the signal peptide/mature polypeptide boundary. The 
stop codons are boxed in red. Figure S3. Fd-GOGAT signal 
sequences. The signal sequences encoded by six plant Fd- 
GOGAT genes are aligned with Clustal V: Wheat (A genome; 
present report and GAJLO 1283868), Brachypodium (BRA- 
DI1G19080), Rice (Os07g46460), Maize (NM_001 1 12223), Ara- 
bidopsis (CP002688), Soybean. (AK245357). The red vertical line 
indicates the exon 1 / exon 2 junction. The blue line indicates the 
end of the signal peptide. Figure S4. Alignment of wheat Fd- 
GOGAT polypeptides. Alignment of the three hexaploid wheat 
Fd-GOGAT polypeptides from the A, B, and D genomes. 
Differences in amino acid sequence are highlighted in yellow. 
Figure S5. Phylogenetic tree of Fd- and NADH-GOGAT 
proteins from diverse species. A selection of available 
GOGAT amino acid sequences from diverse phyla through 
genera were aligned with Clustal W and a phylogenetic tree 
formed by nearest-neighbor analysis. The wheat Fd-GOGAT 
homoeologues are in blue and red for the NADH-GOGAT 
version. Annotated sequences from other organisms have the same 
color coding as wheat. Unannotated sequences are in black. 
Figure S6. Comparison of Fd-GOGAT and NADH- 
GOGAT amino acid sequences. The amino acid sequences 
of the mature Fd-GOGAT from the A-genome is compared to the 
mature NADH-GOGAT from the A-genome by dot blot. 
Matching criteria was 80 match in a 5 amino acid residue 
window. Text SI. Fasta file of the three cv Chinese Spring 
Fd-GOGAT genes. Text S2. Fasta file of the three cv 
Chinese Spring Fd-GOGAT mature proteins. 
(PDF) 

Author Contributions 

Conceived and designed the experiments: DN AB OA AG. Performed the 
experiments: DN OA AG. Analyzed the data: DN AB OA AG. 
Contributed reagents/ materials/ analysis tools: DN AB OA AG. Contrib- 
uted to the writing of the manuscript: DN AB OA AG. 



4. Suzuki A, Knaff DB (2005) Glutamate synthase: structural, mechanistic and 
regulatory properties, and role in the amino acid metabolism. Photosynth Res 
83: 191-217. 

5. Tabuchi M, Abiko T, Yamaya T (2007) Assimilation of ammonium ions and 
reutilization of nitrogen in rice. J Exp Bot 58: 2319—2327. 

6. Canovas FM, .\vila C, Canton ER, Canas RA, Dc La I'orrc E (2007) 
Anmionium assimilation and amino acid metabolism in conifers. J Exp Bot 58: 
2307-2318. 



PLOS ONE I www.plosone.org 



10 



August 2014 I Volume 9 | Issue 8 | e103869 



Fd~GOGAT Genes and Relationship with Grain Protein Content in Wheat 



7. Lccgood RG, Lea PJ, Adcock MD, Hauslcr RE (1995) The regulation and 
control of photorespiration. J Exp Bot 46: 1397—1414. 

8. Ferrario-M6ry S, Hodges M, Hirel B, Foyer CH {2002a) Photorespiration- 
dependent increases in phosphoenolpyruvate cEirboxylase, isocitrate dehydroge- 
nase Eind glutamate dehydrogenase in transformed tobacco pljints deficient in 
ferredoxin-dependent glutamine-alpha-ketoglutarate aminotransferase. Planta 
214: 877-888. 

9. EcrrarioMcry S, Valadicr MH, Godcfroy N, Miallicr D, Hodges M, ct al. 
(2002h) Dhirnal changes in ammonia assimilation in transformed tobacco plants 
expressing ferredoxin-dependent glutamate synthase mRNA in the antisense 
orientation. Plant Sei 163: 59-67. 

10. Laneien M, Martin M, Hsieh M-H, Leustek T, Goodman H, et al. (2002) 
Arahidopsis glt-T mutant defines a role for NADH-GOGA'E in the non- 
photorespiratory ammonium assimilatory pathway. Plant J 29: 347-358. 

11. Sakakibara H, Kawabata S, Takahashi H, Hase T, Sugiyama T (1992) 
Molecular cloning of the family of glutamine synthetase genes from maize: 
expression of genes for glutamine synthetase and ferredoxin-dependent 
glutamate synthase in photosynthetic and non -photo synthetic tissues. Plant Cell 
Physiol 33: 49-58. 

12. Zchnackcr G, Becker TW, Suzuki A, Garrayol E, Gaboche M, et al. (1992) 
Purification and properties of tobacco ferredoxin-dependent glutamate synthase, 
and isolation of corresponding cDNA clones. Eight-inducibilitv and organ- 
specificity of gene transcription and protein expression. Planta 187: 266 274. 

13. Goschigano KT, Melo-Oliveira R, Lim J, Coruzzi GM (1998) Arahidopsis gls 
mutants and distinct Fd-GOGAT genes: implications for photorespiration and 
primary nitrogen assimilation. Plant Cell 10: 741-752. 

14. Avila C, Marquez AJ, Pajuelo P, GanneU ME, Walisgrove RM, et al. (1993) 
Cloning and sequence analysis of a cDNA for barley ferredoxin-dependent 
glutamate synthase and molecular analysis of photorespiratory mutants deficient 
in the enzyme. Planta 189: 475-83. 

15. Boisson M, Mondon K, Tomey V, Nicot N, Laine A-L, et al. (2005) Partial 
sequences ol nitrogen metabolism genes in hexaploid wheat. Theor Appl Genet 
110: 932-940. 

16. Nigro D, Gu Y(^. Huo X, Marcotuli 1, Blanco A, et al. (2013) Structural analysis 
of the wheat genes encoding XADH -dependent glutamine-2-oxoglutarate 
amidotransf erases genes and correlation with grain protein content. PLoS 
ONE 8(9): e73751. doi:10.1371/journal.pone.0073751. 

17. Mdrquez AJ, Avila C, Forde BG, Walisgrove RM (1988) Ferredoxin-glutamate 
synthase from barley leaves: rapid purification and partial characterization. 
Plant Physiol Bioch 26: 645-651. 

18. Bernard SM, BlomMoller AL, Dionisio G, Kichey T, Jahn TP, et al. (2008) 
Gene expression, ccUular localisation and function of glutamine synthetase 
isozymes in wheat iTriticum aestivum L.) Plant Mol Biol 67, 89-105. 

19. (iadaleta A, Nigro D, (iiancaspro A, Blanco A (201 1) The glutamine synthetase 
((jS2) genes in relation to grain protein content of durum wheat. Functional and 
integrative genomics 11: 665-670. 

20. Blanco A. Mangini C>, Ciiancaspro A, Giove S, Golasuonno P, et al. (2012) 
Relationships between grain protein content and grain yield components 
through quantitative trgiit locus analyses in a recombinant inbred line population 
derived from two elite durum wheat cultiv£u-s. Mol Breeding 30, 79-92. 

21. Lawlor DW (2002) Carbon and nitrogen assimilation in relation to yield: 
mechanisms are the key to understanding production systems. J Exp Bot 53: 
789-99. 

22. Triboi E, Triboi-Blondel AM (2002) Productivity and grain or seed composition: 
a new approach to an old problem. Eur J Agron 16: 163-186. 

23. Prasad M, Varshney RK, Kumar A, Balyan HS, SharmaPC, et al. (1999) A 
microsatellitc marker associated with a Q^l'L for grain protein content on 
chromosome arm 2DL of liread wheat. Theor Appl (ienet 99: 341-345. 

24. Khan lA, Procunicr JD, Humphreys DG, Tranquilli G, Schlatter AR, et al. 
(2000) Development of PCJR based markers for a high grain protein content 
gene from Triticum turgidum ssp. dicoccoides transferred to bread wheat. Crop 
Sci40: 518-524. 

25. Zanerti S, Winzeler M, Feulillet C, Keller B, Messmer M (200 1) Genetic analysis 
of bread-making quality in wheat and spelt. Plant Breeding 120: 13—19. 

26. Campbell KG, Finney PL, Bergman CJ, Anderson JA, Giroux MJ, et al. (2001) 
Quantitative trait loci associated with milling and baking quality in a softxhard 
wheat cross. Grop Sci 41: 1275-1285. 

27. Groos C, Robert N, Bervas E, Gharmet G (2003) Genetic analysis of grain 
protein content, grain yield and thousand-kernel weight in bread wheat. Theor 
Appl Genet 106: 1032-1040. 



28. Prasad M, Kumar N, Kulwal PL, Rodcr MS, Balyan HS, et al. (2003) CTLL 
analysis for grain protein content using SSR markers and validation studies using 
NlLs in bread wheat. Theor Appl Genet 106: 659-667. 

29. Zhang L, Spiertz JHJ, Zhang S, Li B, Werf WVD (2008) Nitrogen economy in 
relay intercropping systems of wheat and cotton. Plant Soil 303: 55-68. 

30. Raman R, Milgate AW, Imtiaz M, Tan M-K, Raman H, et al. (2009) Molecular 
mapping and physical location of major gene conferring seedling resistance to 
Septoriatritici blotch in wheat. Mol Breeding 24: 153-164. 

31. Suprayogi Y, Pozniak CjJ, Clarke ER, Clarke JM, Knox RE, et al. (2009) 
Identilieation and \alidaiion of quant itati\'c trait loci for grain protein 
concentration in adapted Canadian durum wheat populations. Theor Appl 
Genet 119: 437-448. 

32. Sun Y, WangJ, CrouchJH, Xu Y (2010) Efficiency of selective genotyping for 
genetic analvsis of complex traits and potential applications in crop 
improvement. Mol Breeding 26: 493-51 1. 

33. Kramer T (1979) Environmental and genetic variation for protein content in 
winter wheat {Tritirum aestivum L.) Euphytica 28: 209-218. 

34. Sears ER (1954) The aneuploids of common wheat. University of Missouri, 
College of Agriculture, Agricultural Experiment Station Bulletin 572: 1—58. 

35. Sears ER (1966) NuUisomic-tetrasomic combinations in hexaploid wheat. In: R 
Riley and K. R Lewis eds. Chromosome manipulations and plant genetics. Oliver 
& Boyd, Edinburgh, pp. 29-45. 

36. Halushka MK, Fan JB, Bendey K, Hsie L, Shen N, et al. (1999) Patterns of 
single-nucleotide polymorphisms in candidate genes for blood-pressure homeo- 
stasis. Nature Genet 22: 239-247. 

37. Till BJ, Burtner C, Comai L, Henikoff S (2004a) Mismatch cleavage by single- 
strand specific nucleases. Nucleic Acids Res 32: 2632-2641. 

38. Till BJ, Reynolds SH, Weil C, Springer N, Burtner C, et al. (2004b) Discovery of 
induced point mutations in maize genes by TILLING. BMC Plant Biol 4—12. 

39. GadaletaA, Giancaspro A, Giove SL, Zacheo S, Mangini G, et eJ. (2009) 
Genetic and physical mapping of new EST-derived SSRs on the A and B 
genome chromosomes of wheat. Theor Appl Genet 118: 1015—1025. 

40. Blanco A, Simeone R, Gadaleta A (2006) Detection of Q_TLs for grain protein 
content in durum wheat, llieor Appl (ienet 1 12: 1 195—1204. 

41. Li H, Yc G, Wang J (2007) A modified algorithm for the improvement of 
composite interval mapping. Genetics 175: 361—374. 

42. Habash DZ, Bernard S, SchondelmaierJ, WeyenJ, Quariie SA (2007) A genetic 
study of nitrogen use in hexaploid wheat in relation to N utilisation, development 
and yield. Theor Appl Genet 114: 403-419. 

43. Fontaine J, Ravel C, Pageau K, Heumez E, Dubois F, et al. (2009) A 
quantitative genetic study for elucidating the contribution of glutamine 
synthetase, glutamate dehydrogenase and other nitrogen-related physiological 
traits to the agronomic performance of common wheat. Theor Appl (jcnet 1 19: 
645-662. 

44. BordesJ, Ravel C, Jaubertie JP, Duperrier B, Gardet O, et al. (2011) Use of a 
global wheat core collection for association Emalysis of flour and dough qusdity 
traits. J Cereal Sci 54: 137-147. 

45. Joppa LR, Cantrell RG (1990) Chromosomal location of genes for grain protein 
content of wild tetraploid wheat. Crop Sci 30: 1059-1064. 

46. Trapp GB, Plank DW, Ganrt JS, Vance CP (1999) NADH-glutamate syntiiase 
in alfalfa root nodules. Immunocytochemical localization. Plant Physiol 119: 
829-837. 

47. Tobin AK, Yainaya T (2001) Cellular eompartinentation of ammonium 
assimilation in rice and barley. J Exp Bot 52: 591-604. 

48. Nakayama T, Archibald JM (2012) Evolving a photosynthetic organelle. BMC 
Biology 10: 35. doi:10T186/1741-7007-10-35. 

49. Miflin BJ, Lea PJ (1980) Ammonia assimilation. In: Miflin BJ, ed. The 
biochemistry of plants, Vol. 5. Toronto: Academic Press. 

50. Sharp PJ, Krcis M, Shewry PR, Gale MD (1988) Location of p-amylase 
sequences in wheat and its relatives Theor Appl Genet 75: 286—290. 

51. Sears ER, Sears LMS (1978) The telocentric chomosomes of common wheat. 
Proceedings of the 5th International wheal genetics symposium. New Delhi: 
Indian Society of Genetics and Plant Breeding, pp. 389-407. 

52. Endo TR, Gill BS (1996) The deletion stocks of common wheat. J Hered, 87: 
295-307. 

53. Van Ooijen JW, Voorrips RE (2001) JoinMap version 3.0: software for the 
calculation of genetic linkage maps. Wageningen: Plant Research International. 

54. Kosambi DD (1944) The estimation of map distances from recombination 
values. Ann Eugenic 12: 172-175. 

55. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: 
molecular evolutionary genetics analysis using maximum likelihood, evolution- 
ary distance, and maximum parsimony methods. Mol Biol Evol 28: 2731—2739. 



PLOS ONE I www.plosone.org 



11 



August 2014 I Volume 9 | Issue 8 | e103869 



